Quantum quenches in the Dicke modeh 
statistics of the work done and of other observables 
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We study tiie statistics of tlie work done in a zero temperature quencli of tlie coupling constant 
in the Dicke model describing the interaction between a gas of two level atoms and a single electro- 
magnetic cavity mode. When either the final or the initial coupling constants approach the critical 
coupling Ac that separates the normal and superradiant phases of the system, the probability distri- 
bution of the work done displays singular behavior. The average work tends to diverge as the initial 
coupling parameter is brought closer to the critical value Ac. In contrast, for quenches ending close 
to criticality, the distribution of work has finite moments but displays a sequence of edge singular- 
ities. This contrasting behavior is related to the difference between the processes of compression 
and expansion of a particle subject to a sudden change of its confining potential. We confirm this 
by studying in detail the time dependent statistics of other observables, such as the quadratures of 
the photons and the total occupation of the bosonic modes. 

PACS numbers: 05.70.Ln, 03.65.Yz 



I. INTRODUCTION 

The study of non-equilibrium phenomena in interact- 
ing quantum systems is one of the most challenging prob- 
lems of modern condensed matter physics because sev- 
eral of the conceptual tools developed to describe physi- 
cal systems in equilibrium ( e.g., the partition function, 
mean field theory, the renormalization group) are not 
readily generalized to non-equilibrium conditions. In or- 
der to make some progress in our understanding of the 
non-equilibrium behavior of these systems, it is impor- 
tant to identify simple paradigms of non-equilibrium pro- 
cesses that may be studied both theoretically and exper- 
imentally. Recently, some progress in this direction has 
been made through the realization of non-equilibrium 
experiments with cold atomic gases loaded in optical 
latticesi'^'^ that, to a good degree of accuracy, are well 
described by simple many-body models such as the Bose- 
Hubbard model."* 

The simplest non-equilibrium process among those 
presently under study is the quantum quench: an abrupt 
change in time of one of the system parameters from an 
initial value A to a final one A'. In a closed system, this 
process corresponds to the preparation of the system in 
the ground state |4'(A)) of an initial Hamiltonian H[X], 
which is then allowed to evolve in time according to a 
final Hamiltonian H[X']. This process is particularly in- 
teresting if some qualitative changes in the state of the 
system between A and A' occur. This was the case in 
Ref. 2, where a gas of bosonic atoms was taken abruptly 
across a quantum critical point from the superfluid to 
the Mott insulating region of the phase diagram. The 
observation of intriguing many-body collapse and revival 
cycles of the two phases in the momentum distribution 
function signaled the high degree of many-body coher- 
ence in the dynamics of these systems^* 

Theoretically, processes assimilable to quantum 
quenches have already been studied a few decades ago in 



a series of seminal papers."' More recently, however, the 
experiments discussed above have inspired an impressive 
and rapidly growing activity on this subject.— 'i'^'^'iSiiii 
Apart from the concrete possibility of testing theoret- 
ical results with experiments, the main motivation be- 
hind this interest has been the shift of focus toward a 
broad class of fundamental issues. More specifically, a 
number of recent studies addressed the extension of the 
concept of universality to the out of equilibrium behavior 
of quantum critical systems subject to either quenches at 
or close to criticality^ or to linear sweeps of the control 
parameter across the quantum critical point. Similarly, 
a great deal of activity is devoted to the search for dy- 
namical manifestations of quantum integrability and to 
the study of the relation between non-integrability and 
thermalizationji 

Looking for a simple and fundamental way to char- 
acterize quantum quenches from the point of view of 
non-equilibrium physics, it was recently observed that 
a quantum quench may be considered in the context of 
basic statistical mechanics as a simple thermodynamic 
transformationi^iii It is thus quite natural to charac- 
terize quantum quenches using standard thermodynamic 
variables: the work W done on the systemj^ii^ the entropy 
S produced,**^ and the heat Q generatedi^i Setting our 
attention on the work done, we note that since a quantum 
quench is not a quasistatic transformation, measurements 
of the work will display fluctuations among different re- 
alizations of the same process. Therefore, for a complete 
characterization of the quench it is necessary to specify 
the full probability distribution of the work done P{W), 
which has to satisfy a number of constraints, such as 
the Jarzynski equalities and the Tasaki-Crooks fluctua- 
tion theoremii^ For quantum quenches the characteristic 
function G{u) = J dWe-'^''P{W) of the distribution 
P{W) was found to be related to the Loschmidt ecYio^^^ 



G{u) = (^'(A)l e''-f^W"e-'-f^[^'l" |*(A)), 



(1) 
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a quantity that emerged previously in the study of X- 
ray singularities in metals^i^ dephasing^i^ and quantum 
chaotic behaviorjii The Loschmidt echo is a simple and 
interesting object: up to a Wick rotation it can be seen 
as the analog of a partition function.^ In addition, its 
direct computation for a prototype quantum critical sys- 
tem, the Quantum Ising chain, revealed that for global 
quantum quenches of the transverse field the presence of 
criticality leads to singularities of the moments of P{W) 
as a function of the quench parameters, while for local 
quenches P{W) itself displays an edge singularity at low 
energies 4 

The goal of the present paper is to move one step for- 
ward toward elucidating and eventually establishing the 
connection between the qualitative features of the statis- 
tics of the work P{W) and the generic characteristics of 
a physical system (e.g., its integrability, the presence of 
a critical point in parameter space). In order to do so, 
it is important to obtain benchmark results for P{W) in 
exactly solvable models, understand its main qualitative 
features, and describe their physical origins. With this 
motivation we study the statistics of the work done in 
quantum quenches in the Dicke model )^^'^^ an exactly 
solvable hamiltonian describing a gas of two level atoms 
interacting with a single electromagnetic cavity mode. 
The Dicke model is an infinite dimensional quantum crit- 
ical system displaying as a function of the atom-photon 
coupling a quantum phase transition between a normal 
and a superradiant phasej ^^'^" 

Focusing on quenches of the coupling constant within 
the normal phase, we show that criticality leaves clear 
signatures in the dynamics of the system and on P(W): 
as the initial coupling tends to the critical point the av- 
erage work done on the system tends to diverge, while 
for quenches ending at criticality the probability distri- 
bution displays an interesting sequence of edge singular- 
ities. We will give a simple physical picture explaining 
this difference, ultimately related to the difference be- 
tween the processes of compression and expansion of a 
particle subject to a sudden change of its confining po- 
tential. We further elucidate these findings by computing 
exactly the time-dependent statistics of observables such 
as the quadrature operators of the photon field and the 
occupation of the bosonic modes. The rest of the paper 
is organized as follows: in Sec.|IT]we introduce the model, 
establish the notations, and discuss the statistics of the 
work P{W). In Sec. Illll we study the statistics of the ob- 
servables starting with the quadratures of the cavity field 
followed by the total occupation of the bosons. Finally, 
in Sec.|TV|we give our conclusions. 



II. THE HAMILTONIAN AND THE 
STATISTICS OF THE WORK DONE 

The Dicke modeliS, describes the coupling of N two 
level systems, such as two level atoms, to a single electro- 
magnetic cavity mode. Its hamiltonian under the condi- 



tion of resonant coupling, which simplifies the succeeding 
analysis, is given by 



N 



X 



N 



-a;oa^a + --=^(a^ + a)((j+ + (j, ), (2) 

^ i—l 



where the Pauli matrices di describe the dynamics of the 
two level systems, while a (at) destroys (creates) a pho- 
ton of frequency ujq- The coupling between atoms and 
photons has been rescaled by 1/VN in order to have a 
well-defined thermodynamic limit A'' +oo. Here and 
throughout this work a natural system of units is em- 
ployed in which action is measured in units of h. 

The Dicke model has two phasesp^ii^ a normal and a 
superradiant one that are separated by a quantum critical 
point at Ac = i^o/^. The transition as A exceeds Ac is 
characterized by the breaking of parity symmetry leading 
to the spontaneous generation of an extensive density of 
photons in the system (a^a) (x TV. 

The Hamiltonian Eq. ^ can be easily diagonalized 
exactly/^ Focusing on the normal phase (A < Ac) in the 
thermodynamic limit it is first of all convenient to re- 
group the Pauli matrices into collective spin operators 
J'^ — J2iLi '^i' where q = z,±. Using now the Holstein- 
Primakoff representation in terms of a bosonic mode 6, 
= 6+6 - N/2, J+ = bWN-b^b, J- = ViV - 6t6 6, 
the semiclassical/thermodynamic limit N — > +oo can be 
taken. One obtains 



H = Loaia^a + b^b) + \{a) + a)(6t + 6) 



(3) 



At this point, the diagonalization proceeds by a standard 
Bogoliubov rotation. The final form of the Hamiltonian 



H[\] = cj+c^c+ -t- w-c^c_ + C. 
Here the eigenenergies are given by 



2A 

uj±{\) woW 1 ± — , 

0^0 



(4) 



(5) 



while the eigenmodes c± can be expressed in terms of a 
and 6 as 

, , , a ± 6 . , , . a'l' ± 6''' , , 

c± = cosh(7-±) — — + smh(r±) — — , (6) 



n/2 



V2 ' 



with tanh(r±) = {ijJ± — lO{))/{lo± -f loq). Finally, the con- 
stant is C = [uj+ + Lu^ — ijJq{N + 2)]/2. The vanishing of 
at Ac = ti-'o/2 is a direct consequence of the presence 
of a quantum critical point. Though we will not consider 
the superradiant phase at length here, we note that for 
A > Ac a consistent thermodynamic limit can be taken 
only after the bosonic operators a, 6 are displaced.— 

In the following analysis we will be interested in quan- 
tum quenches of the coupling constant from an initial 
value A to a final one A'. We will focus on the case where 
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P(W) 





W-5E 



FIG. 1: A typical plot of the probability distribution of the 
work P{W) for a quench from A — 0.5 to A' = 0.9 (the 
delta functions have been Lorentz-broadened for clarity and 
we have set ujq = 1). Notice the sequence of principal peaks 
separated by 2u!+{X') ~ 2.75, each followed by a tail of sub- 
peaks with separation 2lj„(A') ~ 0.63. When A' — Ac these 
subpeaks merge and give rise to an edge singularity at each 
principal peak (see Fig. [2} . 



both A, A' < Ac- From the point of view of the statis- 
tics of the work done, quantum quenches from A < Ac 
towards A' > Ac (or vice versa) are not very interesting 
because the generation of a photon density cx N requires 
a work that scales extensively with the number of atoms, 
while fluctuations are expected to scale like 1/VN, i.e., 
to be highly suppressed in the thermodynamic limit. On 
the other hand, we do not expect major changes in the 
main qualitative results of this paper for quenches with 
both A, A' > Ac: the choice of focusing on A, A' < Ac 
has the only purpose of allowing us to obtain closed an- 
alytic results for the statistics of the work and of other 
observables. 

Let us start our analysis by computing and characteriz- 
ing qualitatively the statistics of the work done in a quan- 
tum quench of the coupling constant from A to A' (A, A' < 
Ac). As stated earlier in the introduction, in order to 
study the probability distribution P{W) of the work W 
done in a quantum quench it is convenient to compute 
its characteristic function G{u) = / dWe"'^'' PiW). In 
terms of the initial and final Hamiltonians H[X] and 
i/[A'], the latter is given by the Loschmidt ccho,^ 



G{u) = (OA|e*^[^l"e^*^[^>|0> 



(7) 



where |0a) is the vacuum of the operators c±. For A, A' < 
Ac the operators diagonalizing the final hamiltonian c± 
are related to the initial eigenmodes c± by a Bogoliubov 
rotation 



c± = cosh(^±)c± + sinh(^±)c]t, 
with tanh(e±) = [uj±{\') - uj±(X)]/[uj±(X') - 



(8) 



^±(A)]. 



From these definitions it is evident that we have to 
compute 



G{u) 



iSEu 



(OA|e-*["+(^')^+'^++'^-(^''^-^-l"|OA), (9) 



where SE is the difference in the ground state energies of 
the initial and final hamiltonians. In order to do so one 
has first to express |0a) in terms of the vacuum |0a') of the 
final eigenmodes c±. Using the definition c±|Oa) =0 to- 
gether with the Bogoliubov rotation Eq. ([5]) , one obtains 

the equation cosh(^-j-)c±|OA) = sinh(^-|-)c]t|OA), which 
implies that 



|0A) = ^+K+]5-[e-]|0A 



where 



(10) 



(11) 



are single-mode squeezing operatorsj^ The state |0a) is 
therefore a squeezed vacuum of the modes c±. In this 
representation the Loschmidt echo takes the simple form 



G{u) = (OA'|^tK+]^IK-]^+K+(^^)]^-[e-(")]|OA'), 

(12) 

where ^±{u) — S,±e 2«i^±(a Using standard formulas 
for the overlap of squeezed states we then obtain 



G{u) = e-'*-^"G+(M)G_(u), 

where 

G±{u) - [1 + n± - e-2*"±(^')"n±]"^ 
Here we introduced the parameters 

[a;±(A')-a;±(A)]2 



n± = sinh (^±) 



(13) 



(14) 



(15) 



which physically represent the average occupation of the 
final eigenmodes c± in the initial ground state |0a). 

From these equations one can immediately deduce that 
the distribution P{W) has the form 



+ 00 



p{w) = 'p+mv^{2i) 

k,l=0 

X S[W - 6E ^2kLu+{X') -2luj.{X')], (16) 



where 



V±i2k) = 



1 (^-\\ 


n± 







(17) 



Qualitatively, the distribution P{W) consists of a series 
of principal peaks separated by 2uj+{X'), each followed by 
a tail of subpeaks describing excited — modes that are 
separated by 2[j_(A') (see Fig. [1]). 

The partial amplitudes 'P± control the weight of each 
peak in P(W). The presence of a quantum critical point 
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P(W) 



W-5E 

FIG. 2: The probability distribution of the work P{W) for 
a quench from A = 0.1 to A' = Ac (here we set ujo = 1). 
A sequence of edge singularities described by Eq. (f29|l and 
separated by 2t<;+(Ac) = 2.75 signals the criticality of the 
system in the final state. 



and its effect on P{W) can be elucidated by studying the 
asymptotic behavior of 'P±{k) for large k. Using Stirling's 
formula, z! ~ ^/27rz^~^^^'^e~^ we obtain 



k 



(2fc)! 



1 



k J 22'=(fc!)2 
from which, for fc 2> 1, one gets 

1 e-'=/«± 



Trfc 



V±{2k) 



1 



1 



(18) 

(19) 
(20) 



The scale C± controls the decay of the corresponding am- 
plitude. Notice now that the vanishing of w_(A) (a;_(A')) 
at the critical coupling implies the divergence of ri_ when 
A ^ Ac (A' Ac). Therefore when A ^ Ac we have 



C- 



4tj_(A) 



oc 



Ac 



Ac - A' 



(21) 



and a similar equation with A' <-> A for A' ^ Ac. The 
presence of a quantum phase transition in parameter 
space is marked by the divergence of the scale controlling 
the exponential decay of the partial amplitude associated 
with the — modes, which are the the ones becoming crit- 
ical at the transition. Notice that when either A or A' is 
exactly at the quantum critical point the partial ampli- 
tude V- decays as a power law, as seen in Eq. p^ . 

Despite the fact that for both quantum quenches to- 
ward the quantum critical point and away from the 
quantum critical point C_ diverges, the physics behind 
these two situations is deeply different. Indeed, while for 
quenches starting close to criticality all the cumulants of 



the distribution P{W) diverge, quenches towards criti- 
cality are marked by the presence of a sequence of edge 
singularities close to each principal peak. At the root of 
this difference, explained in detail below, is the fact that 
for quenches toward criticality the divergence of ^_ is ac- 
companied by the vanishing of a;_(A') as A' — > Ac- Hence 
the energy scale f2_ = 2C_aj_(A') controlling the decay 
of the tails of P{W) remains finite 



(22) 



as A' Ac. On the other hand, for quenches starting 
infinitesimally close to criticality (A — > Ac) we have 



Ac 



4tj_(A) Ac -A' 



(23) 



For quenches starting at criticality (A — > Ac), the diver- 
gence of r2_ implies the divergence of the average work 
done on the system. Physically, this is due to the fact 
that when a quench starts exactly at criticality the mode 
c_ corresponds in coordinate space to a completely delo- 
calized free particle. Every quench to A' ^ Ac describes 
the switching of an harmonic potential tending to confine 
the mode in a finite volume. This process is conceptually 
similar to the compression of a gas occupying an infi- 
nite volume into a finite one: on this basis we expect 
the average work done to diverge as A — *■ Ac. In order 
to see this, we may use the characteristic function G(t) 
to extract directly the cumulants Kn of the distribution 
P{W) using the standard cumulant expansion formula 
G{u) = exp[^,^^^(-z'a)"/n!if„]. Expanding Eq. HID 
to first order we obtain 

(W) ^ Ki = 6E + Lu+{X')7i+ + uj^{X')n^. (24) 

The average excess work {6W) — (W) — SE is then 

K(A')-c.+ (A)]2 , [c._(A')-^-(A)p 



(SW) 



4lu+{X) 



4w_(A) 



(25) 



From this expression it is evident that quantum quenches 
are not reversible processes, W{X A') ^ —W{X' A), 
and that for A — s- Ac one has 



{6W) 



-(A')] = 



1 



(26) 



which diverges at criticality as anticipated. A similar 
divergence is observed in the second moment of P{W), 
which can be easily computed to find 

{{6Wf) = 2K(A')]'n+(l + n+) + 2K(A')]'n-(l + n-). 

(27) 

A diametrically opposite situation is obtained for 
quenches starting with X ^ Xc and ending at critical- 
ity. In this case, the initial potential confining the mode 
c_ is removed in the final state. This process is concep- 
tually similar to the expansion of a gas in vacuum and 



5 



is therefore not expected to be characterized by a diver- 
gence of the average work. However, it turns out that 
the resulting distribution P{W) displays a series of edge 
singularities resulting from the fact that stays finite 
even in the limit A' Ac (see Eq. ([22])). The most el- 
egant way to obtain this result is by observing that for 
u ^ l/a;_(A') we can always approximate 



G-{u) ~ (1 + ir2_u)-i/2, 



whose Fourier transform is 



(28) 



(29) 



Using these expressions for A' = Ac one immediately ob- 
tains 

P{W)^^V+{2k)P-{W-5E-V8LJok), (30) 

k=Q 

which has the expected form (see Fig. [2]) . It is interesting 
to notice that edge singularities in the statistics of the 
work have been previously reported in local quenches of 
the transverse field in a Quantum Ising chain. ^ These two 
examples seem to suggest a connection between criticality 
and edge singularities in the statistics of the work done 
in quantum quenches, a topic which deserves a deeper 
study in the future. 



III. STATISTICS OF OTHER OBSERVABLES 

Let us now continue the characterization of quantum 
quenches in the Dicke model by focusing on the statis- 
tics of observables such as the quadrature operators and 
the occupation of the bosonic modes. In contrast with 
the statistics of the work, which is time independent, the 
statistics of generic observables depends on the time t 
elapsed after the quench. More explicitly, let us consider 
a generic observable Q having eigenstates \n) with cor- 
responding eigenvalues q„. If the initial state before the 
quench is |$o), the probability to obtain q as a result of 
the measurement of Q at time t is 

P{q.t) = 5]|(n|e-*^/*|<i>o)|'<5(q- g„), (31) 

n 

where Hf is the final hamiltonian. Hence, the character- 
istic function GQ{u,t) = J e~^'^"P{q,t) dq of the distri- 
bution P{q, t) is given by the expression 



GQ(7.,0 = (<i>o|e-''^(*)"|$o), 



(32) 



where Q{t) 



AH ft 



Qe 



-iHft 



Let us now compute explicitly this characteristic func- 
tion for two important observables for the Dicke model: 
the quadratures given by 



Xia) 



1 



V2 



(33) 



and the total occupation of the bosonic modes 

N ^a^a + b^b, (34) 

which can be interpreted as the order parameter of the 
transition. 

A. Quadrature operators 

Let us start by computing the statistics of the quadra- 
ture operator X{a) for a generic quench starting at A 
and ending at A' with A, A' < Ac. The appropriate char- 
acteristic function is 

Ga {u, t) = (Oa I e'^[^'l*e-^^(")"e-*^[^'l* |0a) . (35) 

Since the final hamiltonian is diagonalized by the modes 
c±, it is convenient to first express the operator X{a) in 
Eq. (j33p in terms of them. Inverting a transformation 
analogous to Eq. ([6]) (with A ^ A') yields 

1 



a — — = [cosh(r+)c+ -|- cosh(r_)c_ 
V 2 



— sinh(r_|_)c^ — sinh(r_)c_] , (36) 

where the bars mean that the corresponding quantities 
are to be evaluated using the final value of the coupling 
parameter A'. Since the c± operators evolve trivially in 
time according to the final Hamiltonian, we find that the 
characteristic function has the simple expression 

G„(w,t) = (OA|e-*^+("'*)"e-'^-("'*)"|OA), (37) 



where 
X±{a,t) = — - — e 



^ (38) 



with 



A±{a) = [e"''" cosh(f±) - e'" sinh(f±)] . 



(39) 



The next step in obtaining a closed form result consists 
of expressing the state |0a) in terms of the vacuum of the 
operators c±. Using the result of the previous section 
(Eq. Uni)), we obtain 

Gc(u,t) - (Oa' I e"*^+("'*)"e-'^- |0a'), (40) 

where X{a,t) = S\S,±)X±{a,t)S{^±). Finally, using 
the formula 5*^(0 cS'(f) — cosh(^) c 4- sinh(^) c\ we ob- 
tain 



(41) 



with 



-iuj± t 



sinh(^±) 



^;(a)e'"±*cosh(C±)]. (42) 



6 



Since the exponentials appearing in the expression 
Eq. for the characteristic function are standard dis- 
placement operators, taking their vacuum expectation 
value gives^i 



Ga(u,t) = e 



-T-[|')+(a,t)P + l')-(a,t)P] 



(43) 



The statistics of the quadrature operators maintains the 
characteristics it has in the initial state and is always 
gaussian. The only scale characterizing the distribution 
is its variance, which is given by 



{5X{af)^\ri+{a,t)\ +\'q-{a,t)\ 



(44) 



Using this expression we may get further insight on 
the difference between quenches that are driven toward 
criticality and those that start near criticality and then 
are driven away from it. Indeed, focusing on the case a — 
that corresponds to the "coordinate" operator X{a — 
0), we have 



ijj±{\) sin[tj±(A')t]" 
4R(A0F 



I !r,,.\2 f cos[c^±(V)^]^ 

=""1 4c.±(A) 

(45) 

Notice now that in the case in which the initial state is 
close to critical, the closer A is to Ac the more delocal- 
ized the mode c_ is initially. This results in the diver- 
gence of the amplitude of the oscillations of ((5X(a)^) as 
l/\/Ac — A when A Ac. On the other hand, when the 
final coupling constant A' approaches criticality, one is 
describing the physics of an initially confined mode c_ 
that is "released" eX t — Q: it is therefore not surpris- 
ing that for large times t the width of the distribution 
increases linearly with time when A' — Ac, that is. 



{6X{a) 



2\ ^, ^o^-(-^) ^2 



(46) 



B. Occupation number 

The statistics of the total occupation of the bosonic 
modes 



(47) 



turns out to encode similar information. Let us compute 
the associated characteristic function 

Gn{u) = (OaI e'^[^'l*e-*^"e-'^[^'l* |0a). (48) 

First of all we express the operator N in terms of the 
modes diagonalizing H[\']. Using Eq. ([6]) we obtain N = 
N+ + where 

A^_l_ = cosh^(f±)clj_c± + sinh^ {f±)c±cl. 

- sinh(f±)cosh(f±)[(4)2 + (c±)2)]. (49) 

Evolving this operator in time and expressing the initial 
state |0a) in terms of the vacuum |0a') of the operators 



c± as in Eq. [TOl we obtain Gn{u) = Gat^i^ (m)Gjv_ (u), 
where 



Gn± (u) = e 



exp 



'^7.(±)i^,(±) 
j=i 



Here 



(4) 



t ^2 



K3i±) = 



(50) 

(51) 
(52) 



and 



71 (±) = 72(±)* = -m{cosh(2f±) sinh(2^±) 

- sinh(2f±) cosh(2^±) cos[2w±(A')t] 

- sinh(2f±)sin[2w±(A')t]}, (53) 
73(±) = -2m{cosh(2f±)cosh(2^±) 

- sinh(2f±)sinh(2C±)cos[2w±(A')t]}. (54) 

In order to compute the matrix elements in Eq. (jSOp . 
we notice that for both -I- and — modes the operators 
Kj form a closed algebra with commutation relations 
[Ki,K2] = -2K3, [Ki.Ks] = -Ki, and [K2,K3] = K^. 
We may then apply a standard operator ordering theo- 
remSi stating that for this algebra of operators the equal- 
ity 



holds, where 



1,2 



271,2 sinh(/3) 



2/3cosh(/3) -73sinh(/3)' 
1 



cosh(/3) - g sinh(/3) 



2 ' 



(55) 

(56) 
(57) 



with = 7I/3 — 7172 • Using Eq. ((55|) together with 
Eq. ([?T|) . we easily obtain 



exp 



^7,(±)/^,(±) 



(r3(±))'. 



(58) 



Some straightforward algebra now shows that in the 
present case 0^ = —v? and hence a direct computation 
of Fa leads us to 



Gn± {u) 



■iji 



where 



5±(0 = 



\J cos(u) -I- ig± it) sin(M) 
w±(A')^ , wow±(Ao) 



(59) 



cjow±(A) t^±(A')^ 
cj±(A) 



sin2[w±(A')t] 



cos^[cj±(A')t]. (60) 
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The only parameters entering the characteristic func- 
tion are g±{t) {g± > 1). Physically they characterize the 
average occupation of the bosonic modes. Indeed, tak- 
ing the first logarithmic derivative of the characteristic 
function leads to 



(N) 



g+{t)+g_{t) 



1. 



(61) 



As for the average work and the quadrature variance 
{6X{a)'^), the behavior of the occupation for quenches 
starting close to criticality and going toward criticality is 
deeply different. When A' Ac we indeed have that for 
large times 



(62) 



while in the second case (A Ac) the amplitude of the 
oscillations diverges as ^JXc/{Xc — A). 

Finally, we give the result for the full distribution of the 
occupation number A'^. By taking the Fourier transform 
of Eq. [511 one may easily obtain 



+ OC 



P{N) = ^ U{M) S{N - 2M), 



(63) 



where the weights N{M) are given by the finite sums 

M 

. ( k 

k 



U{M) = 



(l+.g+)(l+.g_) 



^ (k - \\ (M - fc - i 



M -k 



^ \ 2k / , X 2(M-k) 

.9+ - 1 \ / .9- - 1^ 



(64) 



to a final one A', both in the normal phase (A, A' < Ac), 
we computed exactly the characteristic function of the 
probability distribution of the work, as well as those 
associated with the statistics of the quadratures of the 
cavity modes and of the total occupation of the bosonic 
modes. We found that while criticality always leaves an 
imprint on the dynamics of the system, and particularly 
on the statistics of observables, there is a deep differ- 
ence between quenches starting close to the critical point 
and those going towards it. In the first case the average 
work (as well as the amplitude of the oscillations of the 
variance of the quadratures and of the average occupa- 
tion) diverges as criticality is approached. In contrast, 
for quenches towards the quantum critical point the mo- 
ments of the distribution of the work stay finite, while the 
distribution itself displays a sequence of edge singulari- 
ties. This occurrence is accompanied by a characteristic 
quadratic growth in time of the variance of the "coor- 
dinate" operator associated with the cavity field, and a 
similar quadratic temporal growth of the average number 
of bosonic modes in the system. We developed a sim- 
ple physical picture explaining the origin of these effects: 
for quenches starting close to criticality the divergences 
are caused by a mode initially delocalized in phase space 
subject to a final confining potential that tends to com- 
press it in a finite volume. The situation is opposite for 
quenches toward criticality: the initially localized mode 
is released at i = 0, and left to spread coherently in 
phase space. On the basis of this general qualitative pic- 
ture, which definitely also holds for quenches with both 
A, A' > Ac, we expect our main qualitative findings listed 
above to apply also in this latter case. 
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